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Abstract. We give a surface integral derivation of the leading-order evolution 
equations for the spin and energy of a relativistic body interacting with other bodies in 
the post-Newtonian expansion scheme. The bodies can be arbitrarily shaped and can 
be strongly self-gravitating. The effects of all mass and current multipoles are taken 
into account. As part of the computation one of the 2PN potentials parametrizing the 
metric is obtained. The formulae obtained here for spin and energy evolution coincide 
with those obtained by Damour, Soffel and Xu for the case of weakly self-gravitating 
bodies. By combining an Einstein-Infeld-Hoffman-type surface integral approach with 
multipolar expansions we extend the domain of validity of these evolution equations 
to a wide class of strongly self-gravitating bodies. This paper completes in a self- 
contained way a previous work by Racine and Flanagan on translational equations of 
motion for compact objects. 
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1. Introduction 

In a recent paper by Racine and Flanagan 0, henceforth referred to as paper I, 
explicit post-l-Newtonian equations of motion for the center-of-mass worldline of an 
astrophysical body interacting with other bodies have been derived from a surface 
integral method. These equations of motion are valid for a large class of strongly self- 
gravitating objects and include the effects of all mass and current multipole moments, 
which are defined by surface integrals of the post-Newtonian gravitational potentials. 

The main purpose of this paper is to pursue the development of the program 
initiated in paper I in two ways. First, using a surface integral method pioneered by 
Einstein, Infeld and Hoffmann [2] and later used in many other references [e.g. Landau 
and Lifshitz j^j, Thorne and Hartle we give a self-contained derivation of leading 
order evolution equations for the spin [Eq. lfHUj) below] and energy [Eq. (jll4j) below] of a 
given body. These two evolution equations have been derived by Damour, Soffel and 
Xu 5J for the case of weakly self-gravitating bodies. We prove here their validity for 
compact objects including neutron stars and quiescent black holes. These equations were 
needed in paper I to obtain final expressions for the translational equations of motion 
and thus the present sequel article completes this computation. Second, we construct an 
explicit expansion of the general solution of a certain post-2-Newtonian field equation 
in a region between two concentric spheres, using a standard method of post-Newtonian 
theory jH]. In this paper we use our solution of a particular post-2-Newtonian field in 
surface integral computations, as part of the derivation of the evolution equation for 
the energy. In addition the expansions of post-2-Newtonian fields outlined in this paper 
could be used in principle to give a more self-contained derivation of the equations of 
motion of paper I. In paper I, an argument invoking globally weak fields was used to 
deduce the form of the laws of motion for strong-field sources. Here we show how to 
rigorously derive, in principle, the laws of motion [Eq. (4.3b) of paper I or Eq. (4.21b) of 
[5]] without invoking any results obtained in the context of globally weak fields. These 
computations would follow the general outline of sectional in which the energy evolution 
equation is derived. 

We first give our conventions below in section 11.11 summarize the formalism in 
section |2j detail in section 01 a construction of a general solution to the relevant post-2- 
Newtonian field equation, then derive, in section 01 the evolution equation for the spin 
and end with the derivation of the evolution equation for the energy in section 

1.1. Notation and conventions 

We use geometric units in which G — c — 1. We use the sign conventions of Misner, 
Thorne and Wheeler [S]; in particular we use the metric signature (— , +, +, +). Greek 
indices (/z,z/ etc..) run from to 3 and denote spacetime indices, while Roman 
indices (a, b, i,j, etc..) run from 1 to 3 and denote spatial indices. The spacetime 
coordinates will generically be denoted by (x°,x l ) = (t,x l ). Spatial indices are raised 
and lowered using 5y, and repeated spatial indices are contracted regardless of whether 
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they are covariant or contravariant indices. We denote by n l the unit vector x l /r, where 
r = \x\ = ^JbijX l xK 

When dealing with sequences of spatial indices, we use the multi-index notation 
introduced by Thorne J9j as modified slightly by Damour, Soffel and Xu [TU] . We use L 
to denote the sequence of I indices a\a 2 ■ ■ ■ cli, so that for any /-index tensor§ T we have 

Tl = 7o iaa ,„ 0j . (1) 

If I = 0, it is understood that Tl is a scalar. If I < then Tl = 0. We define L — 1 to 
be the sequence of I — 1 indices aia 2 ■ ■ ■ a;_i, so that 

Tl-i = T aia2 ^ ai _ 1 . (2) 

We also define 

T%l = T^ aia2 ai . (3) 

For symmetric tensors we use the convention that the version of and the "l + V 
version of define the same tensor. For example T and T ai are defined to be exactly 
the same tensor. We define K to be the sequence of k spatial indices bib 2 ■ ■ ■ b k , so that 

Tr = T blb2 ,„ bk . 

Repeated multi-indices are subject to the Einstein summation convention, as in SlTl- 
We also use the notations 

x L = x aia2 - ai =x ai x a2 ...x ai (4) 

and 

9l — d ai a 2 ...ai = d ai d a2 ... d ar (5) 
We use angular brackets to denote the operation of taking the symmetric trace-free part 
of a tensor. Thus for any tensor Tl, we define 

T <L> = STF L (T L ). (6) 

where STF^ means taking the symmetric trace-free projection on the indices L. For 
example, if I = 2, we have 

The general formula for the STF projection of a tensor is 

[1/2] 

T<L> = ^ ] C k^(aia2---^a2k-l0-2k^a 2 k+i---ai)jiji---jkjk^ (^) 
k=0 

where [1/2] is the largest integer less than or equal to 1/2 and Sl is the symmetric part 
of T L . The coefficients c l k are given by 

J = , l)k » (2/ - 2fc - 1)!! 

k 1 ; (/-2fc)!(2/-l)!!(2A;)!!' 1 ; 

We also use the notation 

xl = x <l> . (10) 

§ Here by "tensor" we mean an object which transforms as a tensor under the symmetry group SO (3) 
of the Euclidean spatial metric <5y , not a spacetime tensor. 
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2. Summary of the formalism 

The physical systems considered in this paper are systems composed of a collection of 
arbitrarily structured astrophysical bodies. It is assumed that the bodies are separated 
enough, and their internal dynamics slow enough, that the gravitational field in between 
the bodies is well approximated by a post-Newtonian expansion (see paper I for a 
detailed discussion). This assumption allows a wide class of strongly self-gravitating 
objects like quiescent black holes and neutron stars. It excludes however bodies emitting 
very strong bursts of gravitational waves. This restriction applies to strong gravitational 
waves produced by internal motions. Strong gravitational waves due to orbital motions 
are allowed. A more detailed discussion of the class of objects and systems encompassed 
by our analysis is given in paper I. 

The systems we consider are expected to be strong sources of gravitational radiation, 
potentially detectable by modern gravitational wave observatories. For example in 
coalescing binary systems, the waveforms of the emitted gravitational radiation are 
expected to carry lots of information, and full exploitation of the expected observations 
will require accurate theoretical models of the waveforms. One of the main motivation of 
the work reported in paper I and in the present sequel paper is to improve the precision 
of theoretical models of gravitational waveforms emitted by systems of compact objects. 

A formalism to describe such systems to post-l-Newtonian order in orbital velocities 
has been developed in detail in paper I. This section summarizes the important features 
of that formalism and serves as a simplified reading guide for sections II and III of paper 
I. The work reported in paper I is, in part, a generalization of the formalism devised by 
Damour, Soffel and Xu [SJ [TUl HU H2] and Brumbeg and Kopeikin [131 EH US] to strongly 
self-gravitating objects. One of the main ingredients of the framework presented in paper 
I is, as in the DSX work, a collection of coordinate systems with a set of post-Newtonian 
field equations and solutions associated with each coordinate system. There is a global 
coordinate system that covers the region in between and far from all bodies. In this 
global coordinate system, each object is treated on an equal footing. There is also one 
coordinate system per body that covers a shell-like region enclosing that body. Each 
of these coordinate systems is adapted to the enclosed body. The precise mathematical 
meaning of the word " adapted" will be spelled out below. The region of overlap between 
an adapted coordinate system and the global coordinate system is called the buffer 
region of the body. All adapted coordinate systems and the global coordinate system 
are conformally cartesian, satisfy harmonic gauge conditions and cover vacuum regions 
in which the post-Newtonian expansion is assumed to be valid. In this paper however, 
we will only need the gravitational field expressed in adapted coordinate systems. All 
the computations requiring the global coordinate system needed for this paper have 
already been performed in paper I. 

Another important part of the formalism is the definition of a set of mass and 
current mulitpole moments characterizing each body and a set of gravitoelectric and 
gravitomagnetic tidal moments felt by each body. These moments are defined as 
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parameters of potentials that are general solutions to vacuum post-Newtonian field 
equations in a particular gauge. We refer the reader to Appendix E of paper I for the 
explicit definitions of the moments in terms of surface integrals of the post-Newtonian 
gravitational potentials. It is important to note that these moments are invariant under 
a post-l-Newtonian transformation of the time coordinate [cf. Eq.(2.30) in paper I], 
which is a gauge freedom analogous to that of electromagnetism. In the DSX papers, 
mass and current multipole moments are instead defined in terms of compact support 
integrals over the matter distribution. As shown in paper I, our field-based definitions 
of mass and current multipole moments coincide with DSX's multipole moments in the 
case of weakly self-gravitating bodies. The field-based definition is more general since 
it can also be applied to strongly self-gravitating bodies. In that case however, the 
moments cannot be expressed in terms of compact support integrals over the matter 
distribution. 

The respective expansions of the metric g^ v and of the Gothic metric g al3 = yfgg alS 
to the required order || for the computations detailed in this paper are, in the coordinate 
system (sajVa) adapted to body A [cf. Eqs.(4.11) and (4.12) of paper I] 



£00 
9oi 
9ij 



1 



2$ y 



2e 2 [(^ A ) 2 + ^ A ]+0(e 4 



e\ A + 0{e* 
5 l3 U - 2 - A 



1 - 2e^ A + e\ - 2^ + 2($^ + xt) ~ e'xti + 0{e 



(11) 
(12) 
(13) 



and 



s 00 = -£ + 4e 3 $ A - e 5 [8($ A ) 2 - A^ A + xt] + 0(e 6 ) (14) 
q o 1 = £ 3 C a + ^ (15) 

Q t] = -^+e 3 xt,+0(e 4 ). (16) 

e 

The terms of order e A in g^ are not required for the computation of the spin evolution 
equation but they will be needed later for the computation of the evolution equation 
for the energy. In Eqs.([li p -([13| h the parameter e is the standard post-Newtonian 
bookkeeping parameter. Our generic time coordinate sa differs from usual time 
coordinates t by a factor of e: 

i= s -f. (17) 

By a usual time coordinate t we mean a time coordinate with the property that g& — ► — 1 
as e — > 0. The motivation behind this choice of a rescaled generic time coordinate is 
explained in section IIA of paper I. The functions appearing in Eqs.lflTj) - (|T6|) satisfy, 
in vacuum, the following field equations 

V 2 $ A = 0, (18) 
|| See the discussion at the beginning of section IV C in paper I. 
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(19) 
(20) 

(21) 
(22) 

(23) 
(24) 

The general solution to the set of equations (|18 p -([20 )) in a region of the form < r_ < 
\va\ < r + can be written as 



V dVAdy A 



These fields are subject to the following (harmonic) gauge conditions 

d$ A d(f n 
4- h t^V = 



and 



ds A dy\ 



CO 



A' S A, 



n S/r A' 1 



=0 
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i=0 



f-1 



\H-1 



/! 



pn M A (s A ) + 



(2/ + r 



d, 



\Va\ 



(25) 
(26) 



(l + l)(2l + 3f L[SA) 

^^f^ n M A (s A )d L 1 -^ - I - *>f (*a)] - ^^f^y "^Mltf [ , (27) 

where 



Hi = ZfjL (28) 

-A 



^ = ^< A l>- (29) 

All tensors, or moments, parametrizing the solutions (j'25|1 - (j27j) are STF on their multi- 
index L only. The moments Z A L and Y^ can be written in terms of fully STF tensors 
by using the following reduction formulas [H HHj 



4 _. A 4/ „„ 2/-1 



~ l _|_ 1 ~ Y+J € ji<ai^L-l>j + 2/ _|_ ^ ^< a i^£-l>' (^0) 

j „a , 2/- 1 

I + l £iko ' 2/ + l' 



A ^ A 2/ — 1 . 

= ^ + T—-r e ^< a i H L-l>j ~ 4 o; , 1 nG L-l> ( 31 ) 
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where overdots mean time derivative with respect to the time argument (here sa)- The 
moments S A and H A introduced above are STF on all their indices. 

The moments n M A and pn M A are respectively the Newtonian and post-Newtonian 
mass multipoles of body A; n G A and pn G A are the Newtonian and post-Newtonian 
gravitoelectric tidal moments felt by body A. The moments S A are the current multipole 
moments and H A are the gravitomagnetic tidal moments. The moments n M A , pn M A , 
n G A and pn G A are defined for I > while S A and H A are defined for I > 1. The /jll and 
vl moments are called gauge moments and give information on the coordinate system 
used. 

The general solution to the field equation (J2*T|) given the expansion (J25|) for the 
Newtonian potential has not been derived in the literature, as far as the author knows, 
and is the topic of section El We therefore do not discuss it here. 

Now there are many harmonic coordinates systems in the region r_ < \ua\ < r+ 
in which the metric has the form (|11 )) -([13 )) . This gauge freedom has been characterized 
in detail by DSX (TO], Kopeikin ^3], Klioner and Voinov ^Ej and in paper I. We define 
the adapted coordinate system to be the coordinate system that takes full advantage of 
this gauge freedom to simplify the solutions (J2HJ)-([2ZJ) in a particular way. It has been 
shown in paper I that this gauge freedom makes it always possible to choose coordinates 
such that the following conditions hold: 



n M A = 0, 


(32) 


P »M A = 0, 


(33) 


A*l = V / > 0, 


(34) 


n G A = 0, 


(35) 


pn G A = 0, 


(36) 


u£ = V I > 1, 


(37) 


Hf = 0. 


(38) 



Equations (JH2J)-(JiIBJ) define the adapted coordinate system. By that we mean that 
after conditions (|3*2*jl - (|3*8*j) have been imposed, there is no gauge freedom left over, except 
constant translations and rotations. Note however that the evolution equations we 
obtain in the end are valid in all mass-centered conformally cartesian harmonic gauges 
since the multipole and tidal moments appearing in the evolution equations are invariant 
under the transformation 

s A = s A + e 4 /?, (39) 

where (3 is a harmonic function. Finally, let us recall here a result from paper I that 
will be needed later on. In paper I expressions for the tidal moments n G A , pn G A and 
were derived [respectively Eqs. (5.25b), (5.30b) and (5.27b) of paper I] in terms of 
the center-of-mass worldlines and all the mass and current multipoles of all the other 
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bodies B ^ A part of the system. Here we will need only the expression for for 
/ > 2: 

oo / -t \k BA 

n °t = E E nJ~ {2k + 21 ~ nM ^5Sr for l > 2, (40) 



where 



and 



B^A k=0 ' BA 



r BA = | V - cra z A \ (41) 



n BA = ( ^ Z B _ c^A) i rBA _ ^) 

The center-of-mass worldlines cm z A are defined in section V C of paper I. 
3. Solutions of 2PN field equations 

The derivation of the energy evolution equation in section|S]requires the general solution 
to (J2U- It can be written as 



Xij = Xi, + E 



-1V+ 1 1 1 

_ / n a 

1=0 



t, — Cijidi —Bij L x L 

l\ r l\ 



(43) 



where x% is a particular solution to (|2T|). with C^l and B^ L being STF on L. Above 
we dropped all the labels A and rename (sa,Da) ~^ C^ 2 ^') to abbreviate the notation, 
since no confusion about the coordinate system can arise: we work exclusively in the 
coordinate system adapted to body A. It is very important to note that the quantities 
CijL and Biji are partially determined by the gauge condition (J2IJ). These moments can 
therefore depend on Z^ and Yn as well as on quantities that are intrinsically of post-2- 
Newtonian order. A similar situation has been encountered at post-l-Newtonian order 
[cf. Eqs. (l3(J|) - (|31j) ] where the moments Zn and Yn parametrizing the gravitomagnetic 
potential depend partially on the Newtonian moments n M L and u Gl, this dependence 
being determined by the harmonic gauge condition 

4$ + did = 0. (44) 

The particular solution x?. is constructed by solving the following elliptic equation 

V 2 f L = r p x L } (45) 

One convenient solution is obtained by using Eq.(A21a) and (A37) of Ref. [U], which are 
respectively 



"(2g-l)H (2*-4g + l)ll 
^ (2q)\(l - 2q)\ (21 - 2q + 1 y A °(^---Oa 2 ^ 1 a 2q ria 2q+1 ...a l ), K^) 



x = r 
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First one starts from the ansatz 

[1/2] 



q=0 



hi 



p+ 2q+2 X § . £. .. 



the log terms being needed to invert ()47|) in the special cases A = I and A = 
The constant b is an arbitrary scale. Applying the Laplacian on both sides of 
requiring that fi satisfies f|45j) yields 



(47) 



(48) 

I - 1. 
and 



a q (l,p) 



l\ 



(2/ -4g + 1)!! 



(49) 



(2q)\\(l - 2q)\ (21 - 2q + 1)!! (p + 2q + 2)(p - 2q + 21 + 3) 
b q (l,p) = 0, (50) 

for q ^ g c , g c being the possible value of q for which (p + 2q + 2) or (p — 2q + 21 + 3) 
vanishes. Note that there exist a q c only if p < 0. For g = q c , a qc (l,p) is undetermined 
and is therefore an integration constant that can be absorbed in the scale b of the log 
term given in We thus set a qc (l,p) to zero by convention. For q = q c , b qc (l,p) is 

given by 

l\ (2/-4g c +l)!! 1 



(51) 



(2q c )\\(l - 2q c )\ (21 - 2q c + 1)!! (2p + 21 + 5)' 
To shorten equations in the rest of the paper, we introduce the following abbreviations 

(5x) 9 L = 5(aia2...ba2 q -\a,2 q Xa2 q +i...cn) (52) 

and 



c q (l,p; r) = a q (l,p) + b q (l,p)lo£ 



(53) 



To solve for Xij, we first write the gradient of the Newtonian potential as follows 



M = E 4^4^^ n M L - I n G iL x L 



H r 2l+3 



l\ 



X1 X <L> 



1=0 



where we used 



,<L> 



1 l (2l + l)Ux< 

r y ' (21 + 1) r^ 1 



(54) 



(55) 



We then have 
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^■* = ££ 



oo oo K T 
0C 00 
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We can then write down the solution 



Xij 
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1 ^ ^ 1) kl 
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km 



T 



(?) 



ijKL 2 
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(56) 



(57) 



n(9) 
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c q (A; + Z, — (2/c + 2/ + 2); r)5i <ak n MK~i>5i< ai n M L 

-i> 



-A- 
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A form of the general solution to Eq. (j21|) is thus 
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+ E 



(-1) 
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(59) 



Again the moments C^x and B i3 -x, are partially determined by the gauge condition 
In the case of weakly self-gravitating bodies, the remaining undetermined pieces of the 
moments parametrizing the homogeneous solution can be expressed in terms of the 
stress-energy tensor of the system. 

In section |S] below, we will need the trace of the term Cjj [cf . Eq. ()43|) above] to 
complete the derivation of the evolution equation for the energy of a body. This term 
turns out to be fully determined by the gauge condition (J24*|) . Substituting (|4"3J) into 
yields 



8 Ay. 



1=0 



/! 



CijL^L ~r:Bij L x L 

r Z! 



<i ~ djXij- 



(60) 



From the field equation (|2ip. we see that djX^j is harmonic. We can therefore expand 
it in the usual form 

'(-l) m 



1=0 



/! 



'QiL^L ~~ ~ \Q iL X <L> 

r Z! 



(61) 
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where Qn and Qn are STF tensors on their multi-index L. Substituting this expansion 
in Eq.([BU|) gives 



E 



-CijLdji— — ~jjlBijjL-ix <L x> 
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Comparing both sides term by term yields, for I > 1 



Ci<L> — -j ( Z iL + Q iL ) . 



Hence we get 



kk 



Q 



kk ■ 



(62) 



(63) 



(64) 



since Zjj equals zero from the gauge condition (|34)1 . The problem is thus reduced 
to identifying the terms in the divergence of the particular solution x^j [cf. Eq.(|57jl] 
that are proportional to x k /r 3 . To complete this straightforward albeit slightly tedious 
calculation, the following formulae, derived using the standard "peeling" identity [e.g. 
Eq.(4.25) of Ref. pj], are useful 

(21 + 3), 
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where Tl is STF on all its indices. The result is 
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1 (1 + 2) „ 
/! (21 + 3) 



(65) 
(66) 
(67) 

(68) 



We will now move on to the derivation of the evolution equations for spin and 
energy. The evolution equation for the spin is simple enough that the material of this 
section is not needed for that computation. It will be needed in the case of the energy 
only. 



4. Spin 

The derivation of the evolution equation for the spin Si presented in this section is 
inspired by the works of Thorne and Hartle [3] and Zhang [Tj\. The method is based 
on the Landau-Lifshitz formulation of the equations of general relativity [H] as 



n 



fiav/3 



,ai3 



16tt 



(-g)T^ u + 



(69) 
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The tensor density r W' OLV ^ is given by 

H Wa>P = gtwgatl _ q^qP^ (70) 

where «/3 = ^/=gg a P. In Eq.flHj), is the usual stress -energy tensor and the Landau- 
Lifshitz pseudotensor T^ u is given by 



q-a/3 



16vr 



,A0 ,m ,A0 ,p i i/ yA/ifl , p ,k + i/Apy ,i/0 |P 



(^^0^0^ + / VXa) + g (2^^^ - 9 aP 9 X ") {2g Vf gar - 9 P .9ur) 0^0^ 

Following Thorne and Hartle we first define the quantity 



• (71) 



Jh = ^- ( f {xw ka0 \ a + n jm ) d 2 ^. (72) 

Here we call an "enclosed angular momentum" for convenience, although it does not 
have a geometric, physical meaning. It is introduced merely as a convenient intermediate 
mathematical tool useful in the derivation of the spin evolution equation. It is a 
meaningful, physical angular momentum (for stationary systems) only if one takes the 
surface S to be a sphere at spatial infinity j.8., which we will never do here. In order to 
derive the evolution equation for the spin of a given body, one must choose the surface 
S to be closed and to lie in the buffer region of that body [see Fig.l of paper I], i.e. 
the surface S must be far enough from the compact object so that it lies in the weak, 
post-Newtonian (or near zone) gravitational field and must exclude all other bodies. 
Substituting Eq.([70)l and the expansion ([T4" )) -([TT) J) into Eq. (|72Jl . we arrive at 

4=it-/ Eijfc (x j Ck,m + SkmQ) d 2 Z m + 0(e 4 ). (73) 

The evolution equation for the spin is obtained in two steps. First, we compute the 
enclosed angular momentum explicitly from Eq.(|73|) and then take its time derivative. 
Next we use the following formula jl], 

j^ = - £e tjk xiT kl d 2 ^. (74) 

Note that the derivation of this formula assumes that the volume enclosed by the surface 
E on a given time slice has the topolgy of M 3 . It excludes, for example, eternal black 
holes. By computing the surface integral on the right-hand side of (|74jl explicitly we 
obtain the evolution equation for the spin of the body. 

We start the derivation by rewriting the expression ()27j) for the gravitomagnetic 
potential Q using to get 



1=0 



1(2/ + !)!!„ x< L > l v ^ <L> 



ZiL n; T i Tr^iLX 



l\ (2/ + 1) r 2 ^ 1 /! 



(75) 



Spin and energy evolution equations for a wide class of extended bodies 



13 



Since Z^ is an STF tensor on its L indices, the following relations hold: Znx <L> = 
Zi<L>% L — Ziix L . Similar relations hold for Yn. The gradient of the gravitomagnetic 
potential f7H)l is then given by 



did = E 



1=0 



1 (2Z + 1)!! 
7! (21 + 1) 



I 



1+2 "V 



Ziji-in 



L-l 



21 + 1 



A+2 



Z iL n 



lr 



l-i 



l\ 



-Y iiL -in 



L-l 



• (76) 



We next substitute (J75J) and (J75J) into (f73j). To perform the integrals, we choose the 
surface S to be a coordinate sphere of radius r. We will drop the S subscripts from now 
on since a specific surface has been chosen. We now have 



J 1 



16n 



(2/ + 1 



(2' + D"(/ + |) WI _ (r , WI 



(2/ + 1) r l + 
!! 1 



^. L7l *i _ r l Y jL n 



e 2 r 2 



16tt 



(2/ + 1) r l+1 3 



kL 



r 2 dQ + 0(e 4 



1=0 



(a + W'+^-c.-iyw' 



(2/ +1) r 



i+i 



dQ + 0(e 4 



(77) 



The integrals appearing in (|77jl are evaluated using the following results, taken from 
Thorne jH] 



1 

47T 
1 

4tt 



2L+1 



rffi = 0, 



— <b n 2L dQ 



21 + 1 



"0(iii2---"i2!-ii2i)j 



(78) 
(79) 



and 



1 

47T 



T K S L n K n L n l dQ 



(l + l) 



-T iL S L if k = l + l 



(21 + 3)!! 

= if |jfe-Z| ^ 1. (80) 

Note that the above equation is true only if Tl and Sl are symmetric trace-free tensors. 
We then obtain after some algebra 

J i = e 2 ^e ijk Z kj + 0(e A ). (81) 
Using Eq. ()30|) . we finally obtain 

,r = e 2 S l + 0(e 4 ). (82) 
Equation (J74j) then allows us to write the spin evolution equation as 

Si = ~^j> e tjk x^T km r 2 n m dQ + 0(e 2 ). (83) 

The space-space piece of the Landau-Lifshitz pseudo-tensor that enters Eq.(JEHJ) is 



T km = |- (d k $d m <f> - -S^d&fyQj + 0(e 4 ). 
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The gradient of the Newtonian potential (|25|) can be written as 



1=0 



n 



14 



55) 



where 



V? 



(2/ + !)!!„ 



4+2 



Mr 



and 



„_ (2Z + 1)!! , 
^iz - , +3 M iL + r G iL . 



(86) 



(87) 



Substituting (J85J) into (J%3|) and then putting (J%3j) into (JE3J), we obtain, dropping terms 
of 0(e 2 ), 



Si = -r 3 e 



ijk 



oo oo „ y -i 

p =o g=o ^ L 

+v* P (#| - V*) r*W> 



dVL 



CO CO . „ 

^ E E ^ f Gtf? - 

(?=0 J 

V 1 ^ + y* ( a y*_y*\ 

^q\(q+l)\(2q + 3)\\ VkjQ[qVQ Vq) 



~y ^ijk 



+E 

p=0 



"5™ ^ijk 



oo 



(88) 



The sum over q in the second to last line above does not contribute since V^q is 
symmetric in j and k. Substituting in (}88|) the expressions PH) and (jH3) for V& and 
t^f respectively, we get 



5; 



E 



1 



+ r l 



Z!(2Z + 3)!! 
(2Z + 1)!! 



(2Z + 1)!! 



4+3 



(Z + 2) n M fcL n M iL + (Z + l)r 2L n G jL »G kL 



(Z + l) n M fcL n G, L -r' 



(2Z + 1)!! 



4+3 



(l + 2) n M jL n G kL 



E 



Z!(2Z + 3)!! 



(2Z+1)!! 



4+3 



^ (Z + 2)"M fcL n M jL + (Z + l)r 2Ln G jX n G' fcL 



(2Z+ 1)!! 



4+3 



(2Z + 3) n M jL n G kL 



(89) 
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The first two terms inside the brackets of Eq. (J89|) do not contribute since they are 
symmetric in j and k. We are left with 



s t =J2 r%J ' L'( 2/ + 1,!! 



Z!(2Z + 1)H 

1=0 K ' 



r l K - J " n M T n G 

' U+3 m 3 L U 



r 



kL 



oo . 

= £-e iJ -j I M jL tt G U ; I (90) 

1=0 

which is DSX's evolution equation for the spin of a given body [Eq. (4.21c) of [3]]. Note 
that the I = term of the sum in does not contribute since n Mj = in an adapted 
coordinate system. We can use the expression ()40|) for the tidal moments u Gl in terms 
of the mass multipole moments of the other bodies of the system to get, restoring body 
labels, 

00 00 i i\k „BA 



St = e ijm ^ ^ ^ ^-(2Z + 2k + 1)!! n Mj L "M*^g^. (91) 

B^A k=0 1=1 ' ' BA 

For a system of Newtonian bodies, the evolution equation (j9"Uj) for the spin of a given 
object can be computed directly by using the Newtonian stress-energy conservation 
equation, which is essentially how DSX derived it. For such bodies, the quantity Sf 
reduces to the usual Newtonian total angular momentum of body A [S], i.e. 

Sf = / d 3 x peijkXjVk, (92) 
J A 

where p is the rest energy density and v % is a Cartesian component of the fluid velocity, 
all quantities being expressed in the adapted coordinate system of that body. One can 
therefore say that Eq. fffif extends the validity of the classical precession of equinoxes 
[in the sense that the precession is caused by a coupling of mass multipole moments 
to an external tidal field and that the coupling has the same functional form as in 
Newtonian physics] to a large class of compact objects. It is important to remember 
that (}9"Tj) is a relation between multipole moments defined as parameters of the metric 
and that for generic compact objects, it is not possible to express these moments in 
terms of volume integrals over the matter distribution such as (J9"2|) . Note also that 
the computations of this paper are not carried out to high enough order to reproduce 
the standard post-Newtonian spin-spin and spin-orbit precession effects. A derivation 
of these effects from a surface integral method requires a more accurate definition of 
spin, again as a parameter of the metric, than the one used in this paper. One requires 
the metric up to post-2-Newtonian order to define the spin of the body in order to 
display these effects in our framework. A systematic analysis of this problem requires 
a generalization of the theory of reference frames to 2PN order and is thus beyond the 
scope of this paper. However if such a framework were available, the method used here 
to derive the precession of equinoxes for compact objects should yield the standard spin- 
spin and spin-orbit precession effects, as well as new terms involving mass and current 
mulitpole moments. This will be the topic of a further paper. 
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Lastly let us recall that in paper I, another set of mass multipole moments M.^ has 
been introduced when writing down the explicit equations of motion for the worldline 
of a given body. They are basically mass moments defined in a coordinate system that 
does not rotate with respect to distant stars. However the fractional difference between 
M.^ and n M^ is of order e 2 so n M£ or A4f; can be used interchangeably in (jHH). 



5. Energy 

In this section we derive the evolution equation for the energy (or mass; we use these 
words equivalently in this paper) n M + e 2 pn M of a given body of the system. As 
explained in section III.B of paper I, the Newtonian contribution n M to the mass is a 
constant, so the evolution equation will involve only the post-l-Newtonian contribution 
pn M. It is important to be clear about what is meant here by "Newtonian mass" and 
"post-l-Newtonian mass". These masses are defined to be the monopole terms of the 
homogeneous solutions for the potentials $ and ip respectively [cf. Eqs. (J2SI) and (|27j) ]- 
The name "post-l-Newtonian" mass for the parameter pn M comes from the fact that ip 
first appears in the metric at post-l-Newtonian order. More details are provided below. 

The derivation of the energy evolution equation follows a similar path to the 
derivation of the spin evolution equation presented in the previous section. First, we 
define an enclosed energy as follows [I] 

Ez = ^-<fH 0a \d%. (93) 

167T J s 

As was the case for J^, the quantity is only a mathematical tool. We choose to call 
it "enclosed energy" for convenience. We do not attribute to it a physical meaning. It 
represents a physical energy only if one takes S to be a sphere at spatial or null infinity 
which we will never do here. The conservation law that contains the information on 
the evolution of the enclosed energy is jl] 



£ s = - j T 0j d%. (94) 

Substituting the definition (fTUj) and the expansion (fTlj) - (fll)j) into (J9l!j) . we obtain 

Ex = -^-£ [e 2 Ad^ + e% + e% - 8$ 2 - X kk) }d% + 0( £ 6 ). (95) 

This expression involves the post-2- Newtonian field whose expansion in terms of 
STF moments has been discussed in section 01 We again take S to be a coordinate 
sphere of radius r. This allows us to simplify some terms in the calculation of the 
enclosed energy by using 

^/(/•;/' jr/->; ; = r 2 ^- I f{rn l )d£l, (96) 



where / is any function and n l = x l /r as usual. We then have, dropping S subscripts 
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The first two integrals have similar structure and we will therefore consider them 
together as follows 



h =±j> (e 2 <S> + e 4 ij) dQ. 



(9* 



From Eqs.()78J) and ()79|h we see that only I = terms of the expansions ()25|) and 
can contribute to I\. Since n M = 0, n G = and pn G = 0, the surface integral I\ is 
equal to 

16tt 



h 



(e 2 n M + e 4 pn M) 



The next surface integral in Eq. (j97|) is 

J 2 = - e 4 8 j $ 2 dft. 
The square of the Newtonian potential is 



(99) 



(100) 



* 2 = £ 



oo K T 

n n 



k,l=0 



kill 



(2k -1) 



(21 - IV r l 

, ~*M K «M L + 2—*M K »G L 



r 



k+l+2 



+ n G K "G L r 



k+i 



(101) 



From Eqs. (j78j) and (f?9*|) . only the k = I terms survive the surface integration. The result 
is 



oo 
1=0 



" ; " n M L n M L + — — — n M L n G L + 



,21 



(21 + 1)V 2Z + 2 



(2/+ l)r 



(2Z+1)! 



n/nf n/~* 



(102) 



We next skip to the integral in ()97)1 involving the gravitomagnetic potential 

J 3 = £ 4 j> rijCj dtt. (103) 
Only I = 1 terms of the expansion (|27|) contribute, giving 

= 0, 



(104) 



as both terms vanish by our choice of coordinates [cf. Eqs. (j32| - (j38|) ]. Had they not been 
eliminated by a coordinate choice, these terms would have been canceled by appropriate 
contributions in Ii [coming from the gauge moments \i and v\ that would have appeared 
due to a different choice of coordinates. The last remaining surface integral in Eq. (|9*7j) 



is 



h = ~£ 4 i> Xkk dVL. 



(105) 
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Substituting Eq.([43)l in the above yields 

h = 4tt£ 4 {^C kk + B kk j - e A j xl k dQ. 

From Eq.(|57j). the trace of Xij is 
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(106) 



oo oo [(fc+0/2] 

k=0 1=0 q=0 
oo oo [(fc+0/2] 



^) q KL T {q) 
kill jjKL 

2^(k + l-2q)\ 



„K-QL- 



k=0 1=0 q=0 



(107) 



In the integral of xtk over the unit sphere, the only terms that can contribute must have 
q = k = I, otherwise the integral vanishes due to orthogonality of {n <L> } or the fact 
that Tjf KL is STF on K and L separately We then have 



/°° 2 l 



(0 



(108) 



Substituting (J55j) into (|108|) yields after some algebra 



4ir 



■in 



OO 



(2/ + 1) 



1=0 



/!2(/ + l)(2/ + 3)r 2/ + 4 "~ 3lj ±r±jL 



„2l+2 



1=0 



l\ (2Z + 2)(2Z + 3)!! 



n GjL n Gji- (109) 



Putting fIBJ) and (fTTTT^ into En.(fTUfi|). adding together Eqs.® , IfTTKty and (fTUoT) and 
then substituting the result into Eq.(|9*7jl finally yields 



E = e 



oo 

4 Eh 

Z=0 



7(/ + 2)(2/ + l) 



2(/ + l)(2/ + 3) r 2 '+ 3 
7 



„2/+3 



(110) 



2(1 + l)(2/ + 3) 

2(2/ + 3)!! 

We can now make direct use of (|94|) . Taking a time derivative of Eq. ()110|) gives on one 
hand 



E = e 



oo 

4 Ei 

1=0 



7(/ + 2)(2/ + l)!! 1 
(Z + l)(2Z + 3) r 2 '+ 3 

(2/ + 3) 2 + 7 



n M jL «M jL 



(21 + 3)!! 



2Z + 3 n/ 



2(Z + l)(2Z + 3) 



"CjL + "MjL n Gji) 



+ e 4 pn M. 



(Ill) 
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From Eqs. dUD-CED and CD, we have 



Air 



dp (iK, - djQ + 3$Si$ 



0(e e 



(112) 



The computation of the surface integral in the right-hand side of using ()112|) 
is performed with the same methods as used before. There is nothing new in this 
computation and therefore we will not report the details here. The result is 



T 0i d 2 Z 
(1-2) 
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1=0 



/!(2/ + 3) 



7r 



21+3 



(2/ + 1)!! 



(I + 5) n 



, 7(Z + 2)(2Z + 1)!! ■ 



(Z + l) '" J " (Z + l)r 2 '+ 3 

The conservation law (|9*%|) then yields 



(1 



(113) 



pn M A 



oo 1 

=E^ 

(i-2) 



(I + 5) 



(2/ + 3) 2 + 7 

(/ + l)(2/ + 3) 2(/ + l)(2/ + 3) 

< a+3 > 2 + 7 



(Z + l)(2Z + 3) 2(/ + l)(2/ + 3) 

oo 1 



(1 + 1) 



(114) 



which is equivalent to Eq. (4.21a) of [H| in an adapted coordinate system. As we have done 
with the spin evolution equation, we can rewrite Eq. (jll4|) in terms of mass multipole 
moments alone. This form of the evolution equation is, with body labels A restored, 



oo oo 



pn M' 



E EE +»+i>» 

B^A k=0 1=1 



jL M K fc+J+2 
' BA 



(i + 2) n 



, M A a] ( f B n <jKL> 
jL M K k+ i + 2 
I i 



BA 



,BA 



^4(2A; + 21 + 3) n M A L M^r'/""' 1 ' 
[l + l) 



BA V" 1 V r BA 

To gain some physical insight into the meaning of Eq. (|114j) , it is useful to look at the 
volume integral expression for the mass monopole of, say, body A, i.e. n M A + e 2pn M A , 
in the weakly self-gravitating limit. This is [5] 



k+l+3 



(115) 



n M A + e 2 pn M A = [ d 3 x ( T 00 + e 2 pi T 00 ) (l + e 2 x ■ V$) , 

J A 



(116) 



where "T 00 and prf r 00 are the Newtonian and post-l-Newtonian pieces of the time-time 
component of the stress-energy tensor of the body in the adapted coordinate system. 
For a perfect fluid, the Newtonian term T 00 is simply the locally measured rest energy 
density p and the post-l-Newtonian term prf r 00 is equal to p(v 2 — 2$) ^H]- The 1PN 
mass thus includes contributions from internal fluid motion and from the gravitational 
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potential. Note that for systems of interacting bodies, it is the full Newtonian potential 
that enters (jll6j) and so the other bodies directly contribute to the integral expression 
for the mass of body A. One can check, using 1PN hydrodynamical equations of motion, 
that for an isolated system the mass ()116|) is conserved 5j. For systems of interacting 
bodies, the tidal terms in the Newtonian potential due to the presence of other bodies are 
entirely responsible for the time evolution of the 1PN mass. Although expression (jll6|) is 
only valid for a weakly self-gravitating body, one should expect physical contributions to 
the 1PN mass of compact objects to include internal fluid motion, gravitational binding 
and the presence of other bodies. However one must rely on numerical computations 
to study quantitatively the dependence of the 1PN mass on the internal structure and 
the environment of the body since no formula like (jll6J) exists for objects with strong 
internal gravity. 

Similarly to the spin case, the computations presented here are not performed at 
high enough order to display spin interactions in the evolution of the mass. One needs 
a definition of the mass of a body accurate to 2PN order and derive the evolution 
equation for this 2PN mass to see such effects. A systematic solution to this problem 
in our framework requires an extension of the theory of reference frames to 2PN order. 

As discussed at the end of section HI the moments M. A can be used interchangeably 
with the moments n M A in the right-hand side of Eq. (jll5p . since the changes produced 
by such a replacement are fractional corrections of order e 2 . Using the definition of M. A 
given in Eq.(5.39) of paper I, the quantity n M A + e 2 pn M A is equal to M. A . An equivalent 
form for the evolution of the energy, setting now e = 1, is thus^f 



M J 



B^A k=0 1=1 
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km 



-(2Z + 2Jfe+l)! 



Mf L M\ 



B n <jKL> 
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B „,BA n <jmKL> 
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117) 



6. Conclusion 



In this paper, we derived, using a surface integral method, leading order evolution 
equations for the spin Sf and the energy n M A + e 2 pn M A of an astrophysical body part 
of a multibody system. They coincide with the results of DSX in the limit of weakly self- 
gravitating bodies. The derivation is valid for a wide class of compact objects and takes 
into account all mass and current multipole moments. These evolution equations had 
been merely assumed in paper I to obtain the final form of the translational equations 
of motion. This paper completes the results of paper I. As part of the computational 
method we have constructed an explicit expansion of the general solution to a post-2- 
Newtonian vacuum field equation in a region bounded by two spherical shells. This type 
of expansion of post-2-Newtonian fields can in principle be used derive the equations 



% The right-hand side of this equation was denoted T 



t en. 



'zf, c ™zf,Mf,Mi 



in Eq.(1.7b) of paper I. 
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of motion given in paper I without invoking results derived for globally weak fields, 
following the lines of the computation presented in section Future work using our 
expansions of post-2-Newtonian fields could involve, for example, an extension to post- 
2-Newtonian order of the theory of reference frames as presented in paper I, studies 
of gauge invariance of tidal heating and conserved quantities at post-2-Newtonian 
order. A conserved energy to post-2-Newtonian order in terms of the center-of-mass 
worldlines and post-l-Newtonian mass and current mulitpole moments could prove 
useful in astrophysical applications like the gravitomagnetic tidal excitations of r-modes 
in coalescing neutron star binaries [T5] . 
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